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ABSTRACT 



Helioseismic inversions reveal a major discrepancy in sound speed between the 
Sun and the standard solar model just below the base of solar convection zone. 
We demonstrate that this discrepancy is caused by the inherent shortcomings 
of the local mixing-length theory adopted in the standard solar model. Using 
a self-consistent nonlocal convection theory, we construct an envelope model of 
the Sun for sound-speed inversion. Our solar model has a very smooth transition 
from the convective envelope to the radiative interior; and the convective energy 
flux changes sign crossing the boundaries of the convection zone. It shows evident 
improvement over the standard solar model, with a significant reduction in the 
discrepancy in sound speed between the Sun and local convection models. 

Subject headings: convection - Sun: helioseismology - Sun: interior 



1. INTRODUCTION 



With the numerous oscillation modes that have been identified with great accuracy and 
so many well-defined global quantities of the Sun, helioseismology allows us to study the 
internal structure of the Sun in unprecedented detail. The observed acoustic modes depend 
predominantly on the sound-speed profile in the solar interior. By inverting solar oscillation 
frequencies it is shown that the difference of sound speed between the Sun and the standard 



solar model (SSM) is fa irly small, except below the base of the convection zone ( jGough et al. 



19961 : iBasu et al.lll997| ). There is a bump in the sound-speed difference between the Sun and 
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SSM at r ~ 0.67-R. Here r is the distance to the center and R is the surface radius of the 
Sun. This sharp feature is generally attributed to the strong composition gradient caused by 
helium settling. Thus several mechanisms have been introduced to soften the composition 
gradient or other wise smooth the sound-speed trans ition at the base of the convection zone 
(for a review, see IChristensen-Dalsgaard et al.ll201lL and references therein), such as partial 
mixing, early mass loss, and modification of opacity, composition or temperature profile. 
However, all these proposed solutions are somewhat ad hoc; some are not even physically 
realistic. 

In fact, a more reasonable explanation of this sound-speed difference may come from 
a more profound aspect of stellar modeling — the convection the ory. In SSM, convection is 
treated with the mixing- length theory (MLT; lBohm-Vitenselll958l ) . which is an oversimplified 
approach to the fully developed turbulence in stars. One serious deficiency of MLT is that it 
is a local theory. Thus the properties of the convective turbulence in any place of the fluid 
field are assumed to depend only on the local structure at that point. This approximation is 
tolerable in the bulk of the convection zone, but it becomes definit ely invalid in the tran sition 
zone between the convective envelope and the radiative interior (IDeng fc Xiongi 120081 ). 



2. THE NONLOCAL STATISTICAL CONVECTION THEORY 

Although with apparent deficiencies, MLT is still the most widely used convection theory 
in stellar modeling. This is mainly because, besides its simplicity, MLT turns out to be very 
successful in most cases. In solar modeling, the near-surface part of the convection zone is 
significantly superadiabatic. The entropy jump between the atmosphere and the adiabatic 
convective region, is a measure of the efficacy of convection, which in MLT is controlled by the 
mixing-length. This in turn determines the solar radius. The mixing length in MLT can be 
adjusted to give the appropriate entropy jump to obtain the correct radius. In the interior of 
the convection zone, the energy transfer by convection is very effective, and the temperature 
gradient is nearly adiabatic, which hardly depends on a specific convection theory. However, 
MLT predicts a sudden stop of convective motions at the boundary defined by buoyancy 
neutrality (i.e., the Schwarzschild criterion). This is not physically realistic since the fluid 
elements will move past the boundary because of inertia. The local nature of MLT makes it 
inadequate for the treatment of overshooting. 



We have developed a statistical theory of nonlocal convection (IXion g|| 1 9 79l . 1 1 9 8 9l ; I D en g et al. 



20061 ). This theory is derived rigorously from the basic hydrodynamic equations to construct 
the dynamic equations of auto- and cross-correlations of turbulent velocity and temperature. 
The nonlocal effects of convection are represented by the third-order moments. Because of 
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the nonlinearity of the hydrodynamic equations, these moment equations form an unclosed 
hierarchy, and thus require a suitable closure approximation before they can be solved. Here 
a gradient-type diffusion approximation, in which the third-order terms are expressed as the 
gradient of the corresponding second-order terms, is adopted. The theory depends on two 
dimensionless coefficients c\ and C2 which are related to turbulent dissipation and diffusion, 
respectively. They can be adjusted to obtain the correct depth of the convection zone and 
an appropriate extent of the convective overshooting zone. Like the free parameter in MLT, 
they cannot be determined within the theory, but have to be constrained by observation. 

By closing the moment equations at the third-order moments, the nonlocal transport 
of turbulent energy and momentum are included in our nonlocal formulation of convection. 
Therefore the convective motions in one place of the fluid field also depend on the properties 
of the fluid in other regions. In the overshooting zone, nonlocal turbulent diffusion serves as 
the driving mechanis m when both buoyan cy breaking and turbulent dissipation decelerate 
the convective flows (IDeng fc Xiongi 120081 ). Besides our gradient-type approximation , the 
moment equations can also be solved u sing other closure schemes (e.g., I Canute] Il992t for a 
detailed discussion, see lGrossmanlll996l and references therein). 



Our nonlocal theory is more capable of dealing with the nonlocal effects of convection, 
though it is more sophisticated and therefore more complicated in applications than MLT 
is. Nevertheless, the theo ry has been successful in solving prob lems of stellar structure, 
evolution, and oscillations (lXionglll986l ; Kiong &: Dengl 120071 120091 ). 



THE SOLAR ENVELOPE MODEL 



We use our nonlocal convection theory to construct a new solar model, in the following 

Mode l NL. Since a commonly used SSM in helioseismology is Model S of lChristensen-Dalsgaard et al 
( 119961 ). Model NL has been calibrated to the same basic properties of Model S, including the 
photospheric radius R = 6.9599 x 10 10 cm and surface luminosity L Q = 3.846 x 10 33 erg s _1 . 
It is calculated with a uniform c hemical composition. To compare with Model S, we use the 



same ch emical composition (i.e., iGrevesse fc Noelslll993l). We also use the OPAL equation 



of state ([Rogers et al.lll9961 ) a nd opacities (jlglesias fc Rogers! Il996l ). For temperature lower 
than 6000 K, opacity tables of iFerguson et al.l (120051 ) are taken as a supplement. 



To avoid extra complications such as nuclear reactions and chemical evolution in the 
core, our calculations are limited to envelope models with the bottom set at 0.3-R. This is 
well below the bottom of the convection zone, and will not compromise our inversion for the 
structure of the convective envelope. 
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Model NL shows major differences from local models near the boundaries of the convec- 
tion zone, especially in the lower transition region between solar convection zone and radia- 
tive interior. Figure [T] shows the behavior of the temperature gradient V = d In T/d In P, T 
being temperature and P pressure, in this region for both Model NL and Model S. It is clear 
that the transition from the convection zone to the radiative interior in Model NL is much 
smoother than that in Model S. In Model S, V changes abruptly from the adiabatic temper- 
ature gradient V a( j to the radiative temperature gradient V ra d crossing the lower boundary 
of the convection zone. However, in Model NL, V is already sub-adiabatic before reaching 
the boundary, thus there is no abrupt change at the boundary. 

Figure [2] shows the fractional radiative energy fluxes F T /F and convective energy fluxes 
F c /F versus fractional radius r/R for both models. In Model S, the convective energy flux 
F c > in the convectively unstable zone, and F c = in the convectively stable zone. In 
Model NL, F c (or the velocity-temperature cor relation) is positive in the convection zone , 
but becomes negative in the overshooting zone ( jXiong fc Dengll2001t iDeng fc Xiongl 120081 ). 
It is clear in Figure [2] that the convective flux F c changes sign crossing the lower boundary 
of the convection zone. In the overshooting zone, F c is negative and the radiative flux F T is 
larger than the total flux of the Sun. As a result, the temperature in the overshooting zone 
beneath the convection zone will increase. Thus the sound speed in this region is higher 
than MLT predicts, which is in accordance with the results of helioseismology. 

In Model NL, the change of sign of the velocity-temperature correlation happens at both 
boundaries of the convection zone. At the upper boundary, this t heore tical prediction has 
been proved by observation. In the classic work of iLeighton et al.l ( 119621 ). they reported the 
reversal of sign of the brightness/temperature- velocity correlation in the sola r atmosphere. 
The same conclusion was drawn by later observations with new technology (ISalucci et al. 
1994]). The properties of the lower overshooting zone can be verified by numerical simulation. 
Two-dimensional and three-dimensional numeric al simulations of compressible convection 
( IHurlburt et al.lll986t iMuthsam fc Zochlinq 119931) show that in the overshooting zone, the 
convective energy flux becomes negative. iKupkal ( 119991 ) confirmed this consistency between 
numerical simulations and the moment equations. 

In local convection models, the lower boundary of the convection zone is defined by the 
Schwarzschild criterion V = V a d- However, the smooth transition of V in the nonlocal con- 



vection model makes it difficult to give such a clear definition. IDeng fc Xiongl ( 120081 ) argued 
that a proper definition of the boundary of the convection zone should be the place where 
the convective energy flux F c (i.e., the correlation of turbulent velocity and temperature) 
changes sign. The physical meaning of this definition is that the convection zone is the zone 
of buoyant force driving for convective motion, while the overshooting zone is the dissipa- 
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tion zone against convective motion. Moreover, using this definition, the local and non-local 
models with the same depth of convection zone will have similar structures. As illustrated 
in Figure 121 the convectively unstable zone is where the convective flux is greater than zero, 
while the overshooting zone is where the convective flux is smaller than zero. By this crite- 
rion, we designed Model NL to have its lower boundary of the convection zone at 0.7123-R, 



which is consistent with l ocal models and helioseismic results (jChristensen-Dalsgaard et al. 



199ll : lBasu & Antialll997f ). 



4. SOUND-SPEED INVERSION 



Helioseismic inversions are based on a linear perturbation analysis of the oscillation 
equations around a reference model. The differences between the frequencies of the Sun and 
the reference model are related to the differences in the structure between them. 



We use the Aarhus adiabatic oscillation package (jChristensen-Dalsgaardl 120081 ) to cal- 
culate the oscillation frequencies for both Model NL and Model S. In order to compare 
with observations, we use data obtai ned from the MDI e xperiment on the Solar and He- 
lioseispheric Observatory spacecraft (jScherrer et al.lll995l ). Figure [3] shows the scaled fre- 
quency differences Q n jSu„j. betwee n the Sun and models, where Q n \ is the inertia ratio 
(jChristensen-Dalsgaard et al.lll996l ) and v n [ is the frequency with radial order n and spheri- 
cal degree /. Q n i8v n i is largely independent of I. This reflects that the frequency differences 
between the Sun and the models are mainly caused by the effects localized near the sur- 
face. However, the scaled differences for Model S (Figure [3](b)) show two distinct branches 
corresponding to modes that respectively do not, and do, penetrate beyond the localized 
difference in sound speed between the Sun and this model (cf. Figure 0J. 

Since Model NL includes only the envelope above the core, it cannot reproduce the 
frequencies of the deeply penetrating low-degree modes. In order to avoid the influence 
of these deeply penetrating modes, we need to make a selection before carrying out the 
sound-speed inversion to make sure that the lower turning points r t of the modes used in 
the inversion are well above the bottom of the envelope model. Our criterion is r t = OAR, 
which leaves 1824 modes from degree / = 10 up to I = 200 for inversion. 



We then employ the subtractive optimally localized averages method (jPijpers fc Thompson 



19921 ) to invert the observed data using Model NL as the reference model. The results are 
shown in Figure H] as filled symbols. For comparison, the results of Model S using the same 
mode set are shown by the open symbols. The well-known bump in sound-speed differences 
near the bottom of the convection zone (r ~ 0.67i?) using Model S is absent in the inver- 
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sion with respect to Model NL. Since the convection formulation is the only ingredient that 
differs between the two models, this progress must be due to our improved formulation of 
convection. 



DISCUSSION AND CONCLUSION 



We have shown that the solar envelope model constructed with our nonlocal convection 
theory has a smooth transition from the convection zone to the radiative interior. The 
velocity-temperature correlation changes sign crossing the boundaries of the convection zone. 
Using such a nonlocal convection solar model, the major discrepancy in sound speed between 
the Sun and SSM can be removed. 

Although other ad hoc methods may also be employed to adjust the sound-speed pro- 
file, we believe that nonlocal convection is a more physically reasonable solution. Before 
diffusion was included in solar modeling, the discrepancy in squared sound speed between 
the Sun and the MLT model was nearly 2% below the convection zone, because the depth 
of th e convection zone in the model was too small (jChristensen-Dalsgaard et al.lll993l ; iBasu 
20101 ) . Inclusion of diffusion increased the depth of the convection zone but led to a strong 
composition gradient in the present SSM. To soften t his composition gradient, partial mix 



ing in the transition zone was then intro duced (e.g., iRichard et al.lll996l ; iBrun et al.lll999 
Christensen-Dalsgaard fc Di Mauroll2007| ) . However, in our nonlocal convection theory, there 
is no abrupt change of physical quantities at the bottom of the convection zone. The non- 
local effects are suitably parameterized, so the solar model has a convection zone of corre ct 
depth and a smooth transition region below it (jXiong fc Dengj|200ll ; iDeng fc Xiongll2008l ). 



Moreover, the connect ive mixing of chemical elements is very efficient in the overshooting 
zone (jDeng fc Xiongi 120081 ) . Here matter is almost fully mixed, therefore there is no abrupt 
change in composition. Such extended mixing of elements may have important influences on 
the chemical evolution of the Sun because it transports fragile elements such as lit hium to 
the hot interior where they are destroyed by nuclear reactions (jXiong fc Deng! 120091 ) . 



We note that modifications to t he tempe r ature profile resulting from overshooting have 
als o been considered by, for example. iRempell (120041; for a helioseismic analysis of the result, 



sec 



Christensen-Dalsgaard et al. 201 ll ) and Zhang fc Li ( 2012 ); the resulting temperature 



gradient profiles were superficially similar to the one obtained here (cf. Figured]). 

Beside solar sound-speed inversion, the nonlocal convection theory also provides a solu- 
tion to some other long-standing pro blems in stel lar modeling, such as the semi-convection 
conflicts in massive star evolution (jXiongl 119861 ). lithium depletions in late-type dwarfs 
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( IXiong fc Deng||2009l ). t he theoretical red edg e of the classical instability strip, and Mira-like 
pulsation of red giants (IXiong &: Dengl 120071 ) . Therefore, to get a better understanding of 
stellar structure, evolution, and oscillation, it is preferable to treat convection nonlocally. 
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Fig. 1. — Temperature gradient V vs. fractional radius r/R for both Model NL (continuous 
line) and Model S (dotted line). The dashed lines show the boundaries of the overshooting 
zone in Model NL. The insert shows V in the near-surface region. 
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Fig. 2. — Fractional radiative energy flux F T /F (continuous line), convective energy flux F c /F 
(long dashed line), and turbulent kinetic energy flux F k / F (dash-dotted line) vs. fractional 
radius r/R for Model NL. The dotted lines show the corresponding fractional fluxes of Model 
S; and the dashed lines mark the boundaries of the overshooting region in Model NL. 
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Fig. 3. — Frequency differences between the Sun and solar models, in the sense (Sun)- 
(Model), scaled by the inertia ratio Q nl . Panel (a) shows the scaled frequency differences 
between the Sun and Model NL for the 1824 modes used for inversion, and panel (b) shows 
the scaled differences between the Sun and Model S for the same mode set. 
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Fig. 4. — Relative differences in squared sound-speed c 2 between the Sun and models, in 
the sense (Sun)- (Model). The open symbols are for Model S, and the filled symbols are for 
the present Model NL. The vertical error bars correspond to the standard deviations based 
on the errors in the observed frequencies, whereas the horizontal bars give a measure of the 
localization of the solution. The dashed lines show the boundaries of the overshooting zone 
in Model NL, and the dash-dotted line shows where V = V a d- 




